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We explore dark matter mechanisms that can simultaneously explain the galactic 511 keV gamma 
rays observed by INTEGRAL/SPI, the DAMA/LIBRA annual modulation, and the excess of low- 
recoil dark matter candidates observed by CoGeNT. It requires three nearly degenerate states of 
dark matter in the 4-7 GeV mass range, with splittings respectively of order an MeV and a few 
keV. The top two states have the small mass gap and transitions between them, either exothermic 
or endothermic, can account for direct detections. Decays from one of the top states to the ground 
state produce low-energy positrons in the galaxy whose associated 511 keV gamma rays are seen by 
INTEGRAL. This decay can happen spontaneously, if the excited state is metastable (longer-lived 
than the age of the universe), or it can be triggered by inelastic scattering of the metastable states 
into the shorter-lived ones. We focus on a simple model where the DM is a triplet of an SU(2) 
hidden sector gauge symmetry, broken at the scale of a few GeV, giving masses of order 1 
GeV to the dark gauge bosons, which mix kinetically with the standard model hypercharge. The 
purely decaying scenario can give the observed angular dependence of the 511 keV signal with no 
positron diffusion, while the inelastic scattering mechanism requires transport of the positrons over 
distances ~ 1 kpc before annihilating. We note that an x-ray line of several keV in energy, due to 
single-photon decays involving the top DM states, could provide an additional component to the 
diffuse x-ray background. The model is testable by proposed low-energy fixed target experiments. 

PACS numbers: 98.80.Cq, 98.70.Rc, 95.35.+d, 12.60Cn 



I. INTRODUCTION 
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Annihilation of positrons near the galactic center gives 
rise to a narrow 511 keV gamma ray line that was first 
observed in 1972 [l|, and which has been confirmed by 
numerous experiments since then, most recently by the 
SPI spectrometer aboard the INTEGRAL satellite 0. 
The signal has two distinct components, one associated 
with the central region of the galaxy (bulge) and an- 
other with the disk. There is as yet no consensus as to 
a conventional astrophysical origin for these gamma rays 
, which evidently originate from excess positrons an- 
nihilating nearly at rest. The apparent axial symmetry 
of the bulge component is a point in favor for proposals 
of models of dark matter (DM) that decays or annihi- 
lates into low-energy positrons, since DM should be dis- 
tributed symmetrically near the galactic center. 1 How- 
ever early proposals of this sort were driven toward DM 
candidates that were nearly as light as the electron it- 
self Q, since the injection energy of the positrons can be 
no greater than a few MeV Q (see however Q). Mod- 
els of MeV scale dark matter that couples to e + e~ are 
highly constrained by low energy collider data, and are 
not (in our opinion) theoretically attractive. If dark mat- 
ter is the source of 511 keV gamma rays, one will need to 
verify its properties by direct detection or other comple- 



1 In addition, the east-west asymmetry in the disk component 
claimed by is not confirmed by the more recent analysis of 
0|, using more accumulated data from INTEGRAL. 
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FIG. 1: Spectrum of states for metastable dark matter mod- 
els. Curves with arrows indicate the sequence of transitions 
for the endothermic case, \2 — > X3 — ^ Xi (left) and the 
exothermic one, X3 - ► X.2 — >• Xi (right). 



mentary means to make the explanation of the 511 keV 
signal convincing. 

In this work we propose and revisit scenarios in which 
a long-lived excited state of DM with mass at the 10 GeV 
scale can scatter into a nearby unstable state, whose mass 
differs by only a few keV. The unstable state decays into 
the ground state with the emission of a low-energy e + e~ 
pair. The decay can be relatively fast, but the energy 
is only released after the slow process of inelastic DM- 
DM collisions occurs. (See however the purely decaying 
variant described below.) There are two qualitatively dif- 
ferent ways to realize this, depending upon whether the 
metastable state is the middle one, requiring endothermic 
scattering, or the top one, leading to exothermic. The 
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mass spectra and sequence of transitions are sketched 
in figure [TJ The endothermic version, in the context of 
500-1000 GeV DM, was first presented in [l(| to try to 
simultaneously explain INTEGRAL, PAMELA (Payload 
for Antimatter Matter Exploration and Light-nuclei As- 
trophysics) and ATIC (Advanced Thin Ionization 
Calorimeter) [l2[ excess electron observations. It was 
subsequently discussed with applications to direct DM 
detection in [l3j], for ~ 100 keV values of the small mass 
splitting. The viability of the scenario for INTEGRAL 
was further explored in 14], but only in the heavy DM 
regime. 2 

It would clearly be interesting if the DM mechanism 
for the INTEGRAL observations was somehow tied to di- 
rect detection of the DM [l3|, E3 ■ Our exothermic mech- 
anism is partly motivated by ref. [l8|], which proposed a 
model involving only the states \2 and X3 ( m our no- 
tation), as a means of explaining two indications of di- 
rect detection of dark matter, namel y th e long-standing 
DAMA/LIBRA annual modulation 19J, and the more 
recent observation of excess low-recoil events by the Co- 
GeNT (Coherent Germanium Neutrino Technology) ex- 
periment [20l | . Ref. [l|[ showed that DM with a mass of 
= 4 GeV and mass splitting of a few keV could be con- 
sistent with these observations, using the exothermic nu- 
clear scattering X3^ — ► X2N'. Their observation is that 
the shape of recoil spectrum is sensitive to modulations of 
the local DM velocity when the scattering is exothermic, 
and this can explain the DAM A observations. Addition- 
ally the overall rate for the same parameters is correct 
for explaining the excess CoGeNT events. The idea of 
ref. [l8| is related to the inelastic dark matter proposal 
plj . which however assumed the scatterings to be en- 
dothermic rather than exothermic. (See [22J for another 
discussion of exothermic scatterings.) 

There have been several proposals for DM in the 5- 
10 GeV mass range to explain the DAM A and CoGeNT 
observations [23j]. Most recently, ref. [24| showed that 
elastic DM interactions could simultaneously explain the 
DAMA/LIBRA and CoGeNT observations if the dark 
matter mass is near 7 GeV and the cross section on nu- 
cleons around 2 x 10 -40 cm 2 . We will argue that our 
endothermic scenario is close enough to being elastic, if 
the small mass splitting SM23 is of order a few keV, so 
that the same analysis applies. Even though such a small 
splitting has little effect on DM-nucleus scattering, it is 
important for DM-DM scattering in the galaxy, in the 
present case where the DM is lighter than the nuclei in 
the direct detectors. Getting the observed rate of galac- 
tic positrons limits the maximum mass splitting in this 



case, to values somewhat lower than those that would 
strongly affect the direct detection rates. 

An interesting variant of the above mechanisms is to 
assume that the unstable excited state is so long-lived 
that it still has a relic population in the galaxy, and so 
does not need to be produced by DM collisions. 3 This 
version has more freedom, in that the rate of produc- 
ing positrons (via decays into xi e+e ~) does not depend 
upon the small mass splitting SM23, whereas the rate of 
inelastic scattering X3X3 ^ X.2X2 is rather sensitive to 
SM 23 . 

Our proposals fit nicely within the framework of dark 
matter with a nonabelian gauge symmetry in a hidden 
sector, as suggested by [17|, since such DM automati- 
cally consists of multiple states, and small mass splittings 
are radiatively generated if the gauge symmetry is spon- 
taneously broken. The simplest example that contains 
three DM states is a hidden SU(2) gauge sector, where 
the DM is in the triplet representation. After the hidden 
SU(2) breaks, two colors of the dark gauge boson must 
acquire small kinetic mixing q ~ 10~ 3 — 10 -6 with the 
standard model hypercharge Y, 

£m ix = 5>i?rv (1) 

i 

while the remaining one must have negligible mixing to 
keep the long-lived state stable against decays to xi e+e ~ ■ 
The mixing Bi's couple weakly to charged Standard 
Model (SM) particles, and mediate the scatterings with 
nucleons or decays into e + e - . 4 We find that the hidden 
gauge symmetry should break at the 10 GeV scale (re- 
sulting in gauge boson masses of order several hundred 
MeV) to give the right cross sections for DM scattering 
in the galactic center and in detectors. An attractive 
feature of these hidden sector gauge boson masses and 
couplings is that they are in the right range to be di- 
rectly probed by new proposed fixed-target experiments 

We will present our main results first, in section [TTJ 
The remaining parts of the paper supply the many details 
leading to these results. Section IIIII specifies the hidden 
sector SU(2) particle physics models we consider. The 
gauge coupling a g of this SU(2) is calculated in section 
IIVI by the requirement of getting the right thermal relic 
density for the DM. There we also work out the crucial 
relative abundances of the excited states. Section [V] de- 
scribes how the rate and angular distribution of 511 keV 
gamma rays are computed. Here we also summarize what 
is believed about the location of gaseous media in the 
galactic bulge where positron annihilation is supposed to 
take place, in respose to criticisms of DM interpretations 



2 The original idea of excited dark matter assumed that only 
the ground state was significantly populated, so that excitation 
through the > MeV mass gap must occur in galactic inelas- 
tic collisions [l5| . However more detailed computations showed 
that the collision rate is not hig h enough with such a large energy 
barrier to overcome fiol . 1 14l . |l6(l . 



s Earlier work on deca ying DM as the source of 511 keV gamma 

rays can be found in [la . l2. r | . 
4 Except in equation we will rescale the to be the mixing 

parameter of Bi with the photon for notational convenience. 
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FIG. 2: Solid curves: contours of log R e + /Robs (the positron production rate) for exothermic dark matter, in the plane of 
the average DM mass M x and mass splitting 5M23- Heavy contour labeled "0" matches the observations. Shaded regions are 
allowed by DAMA/LIBRA, from ref. Each plot takes a different value of gauge boson mass with ^1,3 given by eqs. 

(|4I5[) . DM halo parameters are given by eq. (0. 



of the INTEGRAL observations in ref. Q. In section EJ 
we explain how the gauge kinetic mixing parameter e\ 
is constrained to get the desired rates for DAMA. Var- 
ious astrophysical constraints are addressed in section 
IVIII Our predictions for the masses and couplings of 
the gauge boson B\ that mediates the DAMA and Co- 
GeNT reactions, relevant for direct laboratory searches, 
are presented in section I Villi We conclude in IIXI The 
appendices give further details about the kinetic equilib- 
rium of the DM with the SM, and the cross section for 
DM annihilation. 



II. MAIN RESULTS 
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where a g = g 2 /An and ft is the average mass of the gauge 
bosons. The value of the constant c g depends upon the 
number of hidden sector Higgs bosons that can be present 
in the final state of \\ — > HH annihilations; it can lie in 
the range c g = (1.7 — 2.5) x for the scenarios we con- 
sider. In the following, we assume the dark Higgs bosons 
are heavier than the DM, which yields the top value in 
this range, hence a larger rate of positron production. 
The value of c g is derived in section HVl 



In this section we summarize our main results. The 
details leading up to them will be given in subsequent 
sections. The relevant parameter space is the average 
DM mass M x , the mass splitting 5M23, and the masses 
/ij of the hidden sector gauge bosons Bi that mediate the 
interactions with the standard model. We fix the larger 
mass splitting to be SMyi = 1.1 MeV so that there is 
sufficient phase space for the decay into xi e+e ~ while 
insuring that the decay products are not very relativistic, 
as required by constraints on the injection energy of the 
low-energy positrons |8[. Larger values of JM12 tend to 
suppress the positron rate, and the direct detection rate 
for exothermic DM, because of greater depletion of the 
excited state abundance, but our results are not greatly 
sensitive to the exact value so long as it is less than a few 
MeV. 

Because the nonabelian SU(2) gauge interactions take 
the form 

5Xi $2X3 + cyclic permutations, (2) 

Bi mediates the transition X3 — ► Xi > e t C U C - The strength 
of the gauge coupling g is fixed by the requirement of 
getting the observed relic density of DM from thermal 



A. Exothermic dark matter 

We first consider the exothermic case where \3 is the 
stable excited state. In figure [2] we plot contours of 
logi? e +/i? bs, the predicted rate of positron production 
at the galactic center versus the measured rate, in the 
M X -5M23 plane, where SM23 is the small mass splitting 
between DM states 2 and 3. The contours are super- 
imposed upon the DAMA/LIBRA allowed region of ref. 
[18l ]. To illustrate the dependence on the gauge boson 
masses, each graph has a different value of fj,2, the mass 
of i?2- In the exothermic case, B2 is the color that has 
negigible mixing with the SM (to avoid \3 ~^ Xi decays), 
and so fi2 does not directly affect the rates of either direct 
detection nor galactic positron production. However, the 
class of models we describe in section Mil predicts rela- 
tions between 1x2 and the other gauge boson masses, 

Mi £ \ff4~f4 (4) 
2 

M3 = — 5M 2 3 + M2 (5) 
a g 

The first condition (jj) depends on details of how the 
Higgs mechanism in the hidden sector gives masses to 
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FIG. 3: Left: in the exothermic case, log of Y^/Ytot, abundance of stable excited state X3 relative to the total DM abundance, 
as a function of M x , for several values of the gauge boson mass {12, with /ii,3 fixed as in eq. (|4I5[) and 8M23 = 5 keV. SM12 is 
fixed at 1.1 MeV. Right: same but with varying 5M 12 = 0.5 - 15 MeV, and fixed fi 2 = 1000 MeV. 



H 2 = 500 MeV 
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FIG. 4: As in fig. [2j but for less cuspy Einasto profile with 
a = 0.08, r s = 7.5 kpc, p Q = 0.42 GeV/cm 3 . 



the gauge bosons; we take the inequality to be satu- 
rated, which helps to increase the rate of X3X3 — > X2X2 
scatterings (since Bi is the exchanged boson). The 
second condition ([5]) arises because the mass difference 
M3 — A* 2 determines the radiatively generated splitting 
SM 23 = -^a g (n 3 - /j, 2 ). 

From figure [2] one observes that larger values of fi 2 help 
to achieve a large enough rate of positron production, up 
to some optimal value p, 2 ~ 600 MeV, beyond which the 
rate starts to slowly fall with p, 2 ■ The rise at small p 2 oc- 
curs because increasing jj, 2 inhibits X3X3 ~^ X1X1 down- 
scatterings in the early universe, lessening the depletion 
of the \3 state. Figure[31 left panel, illustrates this more 
directly, where the relic abundance of the excited state Y3 
relative to that of total DM abundance Ytot is plotted as 
a function of M x for several values of p 2 ■ There is a satu- 
ration Y 3 /Y to t — > ~ 1/3 as p 2 approaches the GeV scale, 
for the fiducial value 6M% 2 = 1.1 MeV of the large mass 
splitting. (The right panel of fig. [3] indicates that this 
saturation would occur at higher values of p 2 if 5M\ 2 
is increased. The rate of X3X3 — > X1X1 increases with 
8Mi 2 due to the larger phase space.) Further increase of 



p 2 beyond the optimal point decreases the positron rate, 
because p\ is an increasing function of \x 2l and the rate 
of X3X3 X2X2 transitions goes like /^ 4 . 

It may seem surprising that the rate of positron pro- 
duction is a decreasing function of the mass splitting 
5M 2 3, since the phase space for X3X3 ~~ ► X2X2 increases 
with SM 2 3- However, so does the exchanged momentum 
that appears in the propagator of the virtual gauge bo- 
son, and this has the more important effect of suppressing 
the amplitude; see eq. (1251) . 

Figure [2] shows some overlap between the desired rate 
of positron production and the DAMA allowed region for 
\i 2 >, 200 MeV. For each point in the M X -SM 2 3 plane, we 
have adjusted the value of £1 to obtain the DAMA defec- 
tion rate assumed by ref. [lg. £3 is taken to be <, 10~ 3 ; 
the results shown are insensitive to the exact value. Con- 
cerning £3, an intriguing prediction of our model is that 
each direct detection of the process X3-^ ~~ ► X2N' must 
be accompanied by the subsequent production of e + e~ 
through the decay X2 ~> Xi e+e ~ (whose rate scales as 
e|), so in principle one could look for the positron in co- 
incidence. However, the lifetime for the decay cannot be 
much less than 10 3 s, as we will show in section IVII CI 
Since the speed of DM in the galaxy is of order 10 _3 c, this 
occurs too far from the experiment to detect the e + e~ 
pair. In fact this lifetime is much longer than the age 
of the universe for £3 <, 10~ 8 , leading to an alternative 
possibility for explaining the 511 keV signal via decays 
of primordial X2 , more about which in section III CI 

The rate of positron production through DM excita- 
tions is sensitive to the density profile p of the DM halo; 
it scales like p 2 evaluated near the galactic center. We 
parametrize the shape using the Einasto profile 

P = Pe exp (~ ((r/r s r ~ (r /r s ) Q )) (6) 

A set of values that are often considered to be standard 
are a = 0.17, r s = 20 kpc, p e = 0.3 GeV/cm 3 , r Q = 8.5 
kpc. These values for a and r s are based upon pure 
dark matter iV-body simulations that do not not take 
into account the effects of baryons in the central region 
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FIG. 5: Solid curves: contours of log R e + /Robs (the positron production rate) for endothermic dark matter, analogous to fig. [2] 
for the exothermic case. Einasto halo parameters are given in l[5]). Columns correspond to gauge boson masses p2 — 250, 500 
and 1000 MeV respectively. Shaded regions are 90% and 99% c.l. preferred DM masses for fitting DAMA/CoGeNT data, from 
ref. 0]. 
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FIG. 6: As in fig. but for less cuspy Einasto profile with 
a = 0.12, r s = 12 kpc, p = 0.42 GeV/cm 3 . 



of the galaxy [27J ■ We do not obtain a large enough rate 
of positron production using these numbers. However, 
there is strong evidence that the halo is much more con- 
centrated (cuspy) near the center than these values indi- 
cate, due to the presence of the baryons, which have the 
effect of contracting the density [28j . Table Q] shows the 
profile parameters for six Milky Way-like galaxies from 
the Aquarius simulation, which have been reanalyzed to 
include baryonic contraction (29j . Furthermore it has 
been argued that the local density may be larger than 
the canonical value by a factor of 1.3 — 2 3(J 31 1. We 
find that the exothermic scenario gives acceptable over- 
lap between the INTEGRAL and DAMA-allowed regions 
only if we adopt a DM halo that is very cuspy and has 
a somewhat large density in the solar neighborhood. We 
take the most concentrated example in table HI 

a = 0.065, r s = 5.3 kpc, p Q = 0.42 GeV/cm 3 (7) 

to obtain fig. [5J This could still be considered a 
conservative choice, since ref. [31| argues for p = 
0.43(11) (10) GeV/cm" 3 . With these error estimates, one 



might reasonably consider pq — 0.6 GeV/cm . This al- 
lows for some reduction of the cuspiness of the halo with 
very similar results, to a — 0.08, r s — 7.5 kpc for exam- 
ple. Moreover, we can still achieve reasonable consistency 
using the same cuspy profile while keeping p Q = 0.42 
GeV/cm 3 ; see fig. [H 



B. Endothermic dark matter 

If \2 is the stable state, then the transitions \2 X3 
are endothermic. The energy barrier would tend to re- 
duce the rate of such transitions compared to the exother- 
mic case, but there are other differences that also affect 
the rate. Most importantly, even though eqs. (14151) are 
unchanged, the roles of the gauge bosons B2 and B3 be- 
come interchanged relative to exothermic DM. ^3 now 
controls the rate of X2X2 — > X1X1 downscattering in the 
early universe, hence the relic density of \i- Because /X3 is 
naturally the heaviest of the three gauge boson masses in 
our model (see section lrTl]) , this means that the endother- 
mic scenario leads to a significantly larger rate of galac- 
tic positrons than the corresponding exothermic one. We 
thus adopt a less cuspy halo profile in this case, [1] 

a = 0.08, r s = 8 kpc, p Q = 0.42 GcV/cm 3 (8) 



TABLE I: Characteristics of the density profiles of the haloes 
in the Aquarius galaxy formation simulations of ref. [2!| . Col- 
umn 1 gives the name of each halo. Columns 2-3 list a, r 3 , 
the parameters of the best fitting Einasto model, in the inner 
region of the galaxy. 
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FIG. 8: Allowed regions of ref. 24\ for DM to explain 
DAMA/LIBRA and CoGeNT events. 



Our findings for the 511 keV signal for endothermic DM 
are illustrated in figure [SJ We have the freedom to choose 
even less cuspy profiles if desired, with some accompa- 
nying decrease in the value of SM23, as shown in fig- El 
using the more moderate parameter values a — 0.12 and 
r s = 12 kpc. Fig. [7] shows how the relative abundance of 
the stable state, depends upon the masses M x , 

fi2 and mass splitting 5M\2- For the examples shown, 
8M23 should be 4 keV to match the direct detection 
rate corresponding to fig. HI in the allowed M x region 
that is shaded in figs. [SJ O 

Similarly to the exothermic case, we fix the value of t\ 
to get the desired direct detection rate, while £2, which 
controls the rate of decay X3 — > xi e+e ~ 1 is a f ree param- 
eter. We assumed £2 = 10~ 3 in fig. [SJ In contrast to the 
exothermic case, the results are somewhat sensitive to 
this choice: taking much smaller values of 62 mildly sup- 
presses the rate due to its effect on the relic abundance 
I2, as we will further discuss in section llV CI 

Ideally, the analysis of ref. [24| should be redone for 
our slightly inelastic scattering to see how the overlap of 
the DAMA and CoGeNT allowed regions of fig. M might 
be modified. (For this reason we display our results in 
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FIG. 9: Data points show angular distribution from galactic 
center of observed 511 keV signal, reproduced from ref. 
labeled curves are predictions from DM scattering (lowest, 
dotted curve) and from decay (upper curves), assuming no 
propagation of positrons before they annihilate. 



the same range of DM masses as in fig. [SJ) We leave such 
an investigation to future work. 



C. Decaying DM scenario 

If the gauge mixing parameter €2 or £3 is sufficiently 
small, then the excited state \2 or ^3 (depending upon 
whether the DM is exothermic or endothermic) which we 
have referred to as "unstable" can be as long lived as the 
universe. Let us denote the stable and "unstable" ex- 
cited states by \s and x«- Instead of being produced in 
XsXs XuXu scattering, \u can have a significant relic 
density and produce e + e~ from its slow decays. In sec- 
tion [VB] w e- show that the correct lifetime for producing 
the observed rate of positrons results from taking 



£2,3 



10" 



(9) 



The exact expression depends upon other parameters and 
is given by eq. (|4ip . In particular, the dependence upon 
the DM halo profile is much weaker for decays than for 
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the inelastic scattering scenarios discussed above. We are 
no longer constrained to consider profiles such as (J7J. 

It is intriguing that for reasonable choices of the halo 
profile, the decaying DM scenario can explain the mor- 
phology of the 511 keV signal without the need to invoke 
propagation of the positrons before annihilation. This is 
in contrast to the inelastic scattering mechanisms which 
localize the positron production much closer to the galac- 
tic center, as we next discuss. 



For the decaying DM scenario, it is possible to fit 
the observed angular distribution without any smoothing 
from positron diffusion. Figure [3] shows several examples 
with Einasto parameters a — 0.11 and a = 0.13 that 
pass through all the error bars. These examples are close 
to the ones given in table [H and so could be considered 
realistic in light of baryonic compression of the inner part 
of the DM halo. 



D. Angular profile of 511 keV signal 

If positrons annihilate before propagating, the pre- 
dicted intensity I e + of the 511 keV signal as a func- 
tion of angle is just a reflection of the DM density pro- 
file, 5 through a line-of-sight integral, whose form depends 
upon whether the positrons were created through scat- 
tering or decay: 



I e + (x) 



d:r 



l.o.s. 



\{o~v) jfa, scattering 
decay 



(10) 



A/ v 



The integral is taken along the x direction, where /? SjU is 
the density of the stable or unstable excited DM compo- 
nent Xs,ui proportional to the total density p, and t u is 
the lifetime of Xu- 

For the Einasto profile ([7]) we considered for scatter- 
ings, p 2 is practically a delta function, and so the signal 
would appear to come from a point source, p 2 falls to e~ 8 
of its maximum value at a distance of rs = r s (2a) 1 ^ a . 
Even for the standard profile with a = 0.17 and r s = 20 
kpc, r% is only 35 pc, which subtends an angle of 0.2°. 
However the observed signal subtends at least 8° @; see 
fig. 13 Therefore if the scattering explanation is correct, 
all of the observed width must be due to propagation. 6 
The distance corresponding to 8° is 1.1 kpc, which may 
be astrophysically reasonable, depending upon the struc- 
ture of the galactic magnetic field and the injection en- 
ergy of the positrons. 

Apart from astrophysical mechanisms of positron 
transport [33l l34j |. which strongly depend upon the 
poorly constrained magnetic field of the inner galaxy, our 
model suggests another way in which this widening could 
occur due to the streaming of X3 (X2 in exothermic case) 
before it decays. If the gauge mixing parameter for i?2,3, 
the hidden gauge boson mediating the %3,2 — > Xi e+e ~ 
decay, is sufficiently small, €2,3 ~ 10~ 7 , then X3.2 is so 
long-lived that it will travel approximately 1 kpc before 
decaying. 



5 In the case of inelastic scattering, there is some mild r- 
dependence of (crv) due to the r-dependence of the velocity dis- 
persion, which we neglect here. 

6 Our results differ somewhat from those of ref. [3^. which as- 
sumed a less cuspy halo. 



III. PARTICLE PHYSICS MODELS 

The simplest example of a nonabelian hidden sector 
model consistent with the observations we discuss has a 
dark SU(2) gauge group under which the DM transforms 
as a triplet. The most general form of the Lagrangian 
that we will need is 



C = \Xa{0 ai - M x 6 ab ) Xb - ■faB^B!? (11) 

\yXa£abXb + £ H igg S (A W , E) 



ElAtD B^Y 
A, a D a z nv 



Here A« and E are triplets and a quintuplet respectively 
of the hidden SU(2). Two such triplets are needed in or- 
der to get the kinetic mixing (fTJ) required for the direct 
detection signal and galactic positron production. The 
mixing parameters e$ = (Aj)/Aj arise when the triplets 
acquire VEVs. These mixing parameters lead to a cou- 
pling eeiCOs9\v of the electric current to Bf, in the re- 
mainder of this paper, we rescale £j to remove the Wein- 
berg angle from this coupling. A third triplet is required 
to get the right pattern of mass splittings for exothermic 
DM. The VEV of the quintuplet E ab gives the large ~ 
MeV mass splitting. 

We studied this class of models previously in ref. [35[ . 
It is convenient to take the triplet VEVs to be mutually 
orthogonal (A„ ) = Si a A a . Ref. [35| shows that it is easy 
to construct a potential leading to this pattern. We fur- 
ther restrict the traceless symmetric tensor E a {, to have 
VEVs only on the diagonal, 

(E) = diag(A — B, 2B, —A — B), (12) 

This alignment can be accomplished by suitable small in- 
teractions between E and the triplets. With these VEVs, 
the masses of the gauge bosons are given by 

g 2 (A 2 + A 2 + 2(A + 3B) 2 ) 
g 2 (A 2 1 + A 2 3 + 8A 2 ) 
pi = g 2 (A 2 1 + A 2 2 + 2(A-3B) 2 ) (13) 



Mi 

i4 



The corresponding mass shifts in the Xa states relative 
to the average mass M x are given by 



5Mx = ~-a g (p 2 + p 3 )+y(A- B) 



6M 2 
6M 3 



a g (/i! + p 3 ) + y{2B) 

a g {pi+ p 2 )+y{-A- B) (14) 
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We have introduced the Yukawa coupling contribution 
in order to explain two different scales of mass splittings: 
SM 23 ~keV, and AM 12 ~ SM 13 ~MeV. This can oc- 
cur if the quintuplet VEVs satisfy A — —3B; then the 
Yukawa term only contributes to the large mass splittings 
and not to SM 23 . Let us assume this to be the case; 
we will presently show how it can come about. Then 
the small mass splitting comes entirely from the one- 
loop self-energy contribution from gauge boson exchange, 
SM 23 = SM 3 — SM 2 = ka g (fj, 3 — [i 2 ). The assumed order 
5M 3 > 6M 2 requires that fi 3 > fi 2 , hence A 2 > A3. The 
most economical choice would be to remove A3 from the 
spectrum altogether, which is permissible if the kinetic 
mixing parameter e 3 is allowed to vanish. In fact for 
the endothermic scenario, this is exactly what we want, 
in order to forbid \i ~^ Xi e+e ~ decays, so that \i can 
be stable. 7 For the exothermic case, it is opposite: we 
need to insure the stability of X3 against decays to Xi, 
hence e 2 must be negligible, while e 3 is needed for the 
X2 — > Xi e+e ~ decays. And for both scenarios, e\ must 
be nonzero to enable direct detection via inelastic X2, 3 
scattering on nucleons. The upshot is that we need all 
three triplets for exothermic DM (although only two of 
them should lead to kinetic mixing), but only two, Ai. 2 , 
for endothermic. 

Now we return to the question of why the quintuplet 
VEVs should satisfy the seemingly fine-tuned relation 
.4 = —3B. Interestingly, the desired VEVs can arise 
from the simple renormalizable potential 

V(E) = A s (tr£ 2 -v 2 ) 2 +/idct£ (15) 

which has three degenerate minima at A = ±3-8 and 
A = 0. (In the absence of the triplet VEVs, this would 
leave one of the three gauge bosons massless, breaking 
SU(2)— >U(1)). We assume that it is possible to design 
small interactions with the triplets that align (E) along 
the diagonal, and which might perturb A slightly away 
from — 3.B. For simplicity we take A = —3B in the fol- 
lowing, so that 5M\ 2 = 6yB. 

Considering the smaller mass splitting, if M x ~ 5 
GeV, the gauge coupling is of order a ~ 10 -4 , and 
the difference in gauge boson masses should be of order 
\8M 23 \/a g ~ 100 MeV. This is consistent with triplet 
VEVs at the scale of ~ 30 GeV, since the gauge coupling 
is g ~ 0.035. To get the correct sign for the mass differ- 
ence, M 3 > M 2 only requires that A 2 > A3 (given our 
assumption A = —3B), so that fi 3 > fi 2 . 

IV. RELIC DENSITY 

In ref. [35| a first attempt was made to compute the 
value of the gauge coupling a g corresponding to the ob- 



7 See section I VII Al though, for further discussion of some sub- 
tleties. 



served DM density through thermal freeze-out. In this 
section we correct and refine that result, taking into ac- 
count some subleading effects, including extra annihila- 
tion channels into dark Higgs bosons, dependence of the 
Born cross section on the DM velocity, and Sommerfeld 
enhancement (3(| at the time of freeze-out. Moreover we 
estimate the amount of dilution of the excited states due 
to downscattering in the early universe. 



A. Annihilation cross section; determination of a g 

Ref. [35J derived the annihilation cross section for 
XX BB by separately considering X1X1 BjBj (with 
j = 2,3) and XiXj ~~ * B%Bj (again with j = 2,3), and 
explicitly averaging over the initial state and summing 
over the final state colors. In this paper, we make sev- 
eral improvements to the previous calculation as well as 
minor corrections. Details of our calculation are given in 
appendix [XJ 

First, we now include the process where XX goes to 
two hidden sector Higgs bosons through exchange of a 
virtual B in the s-channel. Furthermore, as pointed out 
in [37[, velocity-dependence (including Sommerfeld en- 
hancement) may make significant contributions to the 
cross section in some models. Therefore, we include cor- 
rections to order v 2 in the tree-level cross section as well 
as the leading contribution from Sommerfeld enhance- 
ment. 

Moreover, in the present application, the dark matter 
is sufficiently light that its annihilation may take place 
in the broken phase of the theory. Therefore, we have 
computed the cross section for XX — * BB taking into ac- 
count the gauge boson masses. The main effect comes 
simply from the reduction in phase space, which is a fac- 
tor of (1 — /if /M 2 ) 1 / 2 in the cross section, where fa is the 
mass of the gauge boson in the final state. Annihilation 
to light Higgs bosons is also modified only by this factor 
(of course, Higgs bosons heavier than the DM are not 
annihilation products). For simplicity we replace fa by 
the average mass p, of the gauge bosons and light Higgs 
bosons in this part of the calculation. 

The full cross section for annihilation into gauge 
bosons, N 3 Higgs triplets, and N$ Higgs quintuplets, in- 
cluding leading velocity dependence, is 

♦(^-TT-Tr)*]^)"™ 

in the center of momentum frame. The factor (1 + 
7ra g /u rc i) incorporates Sommerfeld enhancement neglect- 
ing the masses of the gauge bosons. This neglect is 
valid for large M x , such that the freezeout temperature 
~ M x /20 is above the symmetry breaking scale. For 
smaller M x , this factor is roughly an upper bound on 
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Sommerfeld enhancement (except very close to a reso- 
nance) and furthermore a g is sufficiently small that Som- 
merfeld enhancement is an unimportant correction dur- 
ing freezeout. Conversion to the rest frame of the cosmic 
fluid introduces an additional correction of order V 2 for 
center of momentum velocity V. 

To compute the relic density, one needs the thermal 
average (crv ie \). Using the Maxwell-Boltzmann distribu- 
tion, we find 



{(JV rc l) = — 



12 M 2 



/317 



lOiVf 



irM x 
T 



1 T 

2ttMZ 



5^3 

2 



257V, 



T 

mZ 



M 2 



(17) 

1/2 



Annihilations go out of equilibrium at a temperature 
given by M x /Tf = Xf = ln£ — ilnln£ with £ 
1 .0 x 10 12 (M/TeV) for triplet DM, giving x f = 20.4 
for M x = 5 GeV and x f = 23.3 for M x = 100 GeV. 
We should equate (|18p at the Tf with the cross section 
needed to match current constraints on the DM den- 
sity. This varies mildly with M v as (cru)o — (3.2 — 
0.24 log(M/GeV))xlO- 26 cm 3 /s >8]. However, the lat- 
ter expression assumes the usual particle content of the 
standard model at the time of freezeout, whereas in our 
model there are three additional gauge bosons and ex- 
tra dark Higgs bosons. This increases both the Hubble 
rate and redshifting between freezeout and the present. 
Thereby, the extra particle content decreases (<rv)o by a 
factor of x/l + {6 + 3N 3 + 5iV 5 )/61.75. Details are given 
in appendix[X] For two Higgs triplets and one quintuplet, 
we find that a good approximation is given by 



1.7 x 10- 



(M x /GeV) 



(1 



7.7a„ 



(18) 



For the range of M x we are interested in, a g is so small 
that the Sommerfeld enhancement factor y/1 + 7.7a g can 
be neglected. Further results for other light Higgs states 
are given m 1X201 



B. Relative density of excited state 



the rate of positron production. The direct detection rate 
for XsN — > XuN' similarly scales like (n s /ntot)- 

To compute the dilution of Xs from downscattering, 
we solve the Boltzmann equation for the abundance 
Y s = n s /s, where s is the entropy density. Defining 
z = SMi s /T, it can be cast in the form \3£ 



^ = _A Y 2 -Y 2 (XT 
dz z 2 I s 1 {Y^ 



(19) 



where A is related to the cross section 04, for XsXs — > X1X1 
downscattering by 



(20) 



except (as we shall describe below) the multiplicity fac- 
tors <?* and g* s that appear in A should retain their z- 
dependence (only explicit factors of T get replaced by 
SM ls ). 

We can simplify this by assuming that the abundance 
of the ground state does not change significantly during 
the depletion of Xs, so Y\ is just a constant. Furthermore 
Yj q /Y' 1 eq = e^ 2 to a good approximation if the DM is 
in kinetic equilibrium with the standard model particles 
(we will discuss this caveat below) . Defining the fraction 
f = Y s /Y 1 , (HU) becomes 



— = -— (f- 

dz z 2 



(21) 



where A = XY\. (However the e~ 2z will be modified 
when we take into account kinetic decoupling effects; 
see next subsection.) To explicitly compute XY\, let 
us parametrize the DM ground state density as ri\ — 
(g*s/g*s,o)£,T 3 , where the g* s factors take into account 
the dilution of n\ as a function of temperature due to 
entropy production after freezeout. Then 



A = 



g* s t5M ls M. 



l.G6*/g7 g* s ,o 



(22) 



where M p = 1.22 x 10 19 GeV and £, = 7 x 10~ 10 GeV/M x 
to match the observed DM density. 

The cross section is straightforward to compute, 
since it is similar to e~e~ scattering, with just two dia- 
grams, exchange of a gauge boson in the t and u channels. 
In the low- velocity limit, we obtain 



Let us denote the stable excited state by % s , and the 
unstable one as x„. At the freezeout temperature, all 
three DM states are equally populated, but if the rate of 
downscattering transitions XsXs X1X1 remains larger 
than the Hubble rate at temperatures below the mass 
splitting SMis, the density n s of the excited state gets 
suppressed relative to n\ of the ground state. The rate 
of 

XsXs ~~ ^ XuXu transitions in the galaxy at the present 



epoch scales with (n s /n to t) (f° r n tot — n\- 



We 



must therefore compute this ratio to accurately predict 



(a^v) = Ana 2 



M 2 v t 



{^+M 2 x v 2 ) 



2\2 



(23) 



where Vt — \/2SMi s /M x is the velocity of xi at threshold 
(when the incoming Xs particles are at rest) and is the 
mass of the exchanged gauge boson. For the exothermic 
DM model, where Xs = X3i * — 2, while for the endother- 
mic case where Xs — Xi-, i — 3- Using this constant cross 
section is actually a conservative estimate, as it is near 
the maximum value of the full velocity dependent cross 



10 



section given in equation (l36l) for 8M\ S ~ 1.1 MeV and 
our typical values of M x , fa. 

To solve the Boltzmann equation, we first tried to 
employ the semianalytic technique popularized in Kolb 
and Turner [40j. Namely, one writes / = e~ z + A and 
linearizes the equation in A for the early time behav- 
ior, giving A = z 2 /(2A), while A' = — Az~ 2 A 2 at late 
times. Integrating the latter equation between the z 
of freezeout, Zf, and infinity gives the final abundance 
= Aqo = Zf/X. The trick then is to appropriately 
determine the value of Zf. One does this by assuming 
that 



(24) 



and then finding the value of c for which this procedure 
best reproduces the full numerical solution. However we 
find that this procedure is not sufficiently accurate for 
the regime we are interested in, where Y s /Y± 0.1 rather 
than the exponentially small values of interest for Y\ it- 
self. There is no fixed value of c that accurately gives the 
same as the numerical result as Y a /Y\ ranges between 
0.1 and 1. Therefore we numerically solve the Boltzmann 
equation in all cases. 

To relate Y s /Y\ to the fractional abundance of the 
stable state to the total dark matter population at the 
present time, we must remember that the unstable state 
Xu is a lso kept in equilibrium with xi until a similar 
freezeout temperature (which is the same in the limit 
that /i2 = /U3). Only at much later times > 10 3 s, Xu 
decays to the ground state. The total abundance of dark 
matter is then Y\ +Y U + Y S . The fractional abundance 
of Xs is therefore 



Ytn\ 



Y./Yi. 



l + Y s /Y 1 +Y u /Y 1 



(25) 



where Y SjU /Yi denotes the values at freezeout, from solv- 
ing the Boltzmann equation. We compute Y u jY\ in ex- 
actly the same way as Y s /Y\, The only difference is the 
exchange of ^2 ^ f°r the- gauge boson mass appearing 
in the propagator of the cross section (|2"5)l . An astute 
reader may wonder whether X3X3 — > X1X1 scatterings 
change the ratio Y u /Y s additionally; however, the cross 
section f[23]) is greatly reduced for SM23 ~ 1 — 10 keV, 
and we find like [l8| that this process freezes out at tem- 
peratures well above the mass splitting SM 2 3- 



C. Kinetic equilibrium 

The preceding discussion of the Boltzmann equation 
assumed that the DM is in kinetic equilibrium until the 
freezeout of downscattering. If this is not the case, the 
relic density of Xs will be smaller than estimated there. 
The reason is that the equilibrium density depends upon 
the kinetic temperature Tk and this redshifts with the 
expansion of the universe Tk ~ 1/a 2 , in contrast to the 




FIG. 10: Scattering of \ on charged particle q that keeps DM 
in kinetic equilibrium. 



temperature of particles that are still coupled to the ther- 
mal bath, T ru 1/a. To get some sense of the size of the 
effect, we can follow the analytic procedure for an ap- 
proximate result, even though in the end we solve the 
Boltzmann equation numerically. 

If T d is the kinetic decoupling temperature, then Tk = 
T 2 /T d for T < T d . Let z d = 5M ls /T d . Then the term 
e~ 2z in p^|) must be replaced by exp(— 2 max(z, z 2 /z d )). 
Following the semianalytic approach described above, 
one finds that eq. is replaced by 



A(z/) 



c\z d 



z}/zd 



(26) 



which can be rewritten as Zf — (z d \a{c\z d / z^)) 1 / 2 . As 
a consequence the relic abundance of x.s is suppressed 
by ^fz^ in this case. Thus it is preferable for kinetic 
decoupling to occur after the chemical freezeout of % s , 
for maximizing its relic density. 

The principal interaction for maintaining kinetic equi- 
librium with the SM is the electron scattering diagram 
shown in fig.[T(J] The rate for this process is computed in 
appendix[Bj The decoupling temperature as a function of 
e (the kinetic mixing parameter for whichever gauge bo- 
son is exchanged) is shown in fig. [TT] for the case M x = 5 
GeV and /i = 100 MeV. This can easily be generalized to 
other DM and gauge boson masses by noticing that the 
rate scales like a g e 2 /[i 4 and a g is proportional to M x . 
Hence the scaling of e in fig. [TTJ For lower values of e 
than shown in the figure, the relation extrapolates to a 
power law, 



T 



10 MeV 



1.2 x 10- 



-2/3 



(27) 



In reality there are two transitions with two different 
mass splittings that can maintain kinetic equilibrium, 
since we also have the xi e X3 e reaction with the small 
mass splitting 6 M23. We compute the decoupling tem- 
perature for both reactions and take the smaller of the 
two as the true T d . Roughly speaking, only the larger of 
the two e's is therefore relevant for kinetic equilibrium. 
Figure [TT] shows that there is a weak dependence upon 
SM with the large mass gap giving a bigger effect. There 
is also dependence upon the gauge boson masses. 
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FIG. 11: Decoupling temperature for process of fig. [T0]as a 
function of gauge kinetic mixing parameter, for several values 
of the large mass splitting, and for DM mass M x = 5 GeV 
and gauge boson mass fi = 100 MeV. 



For the exothermic model, the coupling e\ is large 
enough so that e 2 is practically irrelevant for kinetic equi- 
librium. This is illustrated in figure [12] (left panel), which 
shows that contours of log R e + /i? b s hardly change be- 
tween e 2 = 10 _1 and 10~ 10 . (The example shown is 
for fi2 = 250 MeV; for larger the dependence is even 
weaker.) For the endothermic model, e\ is smaller and 
so £3 can have a bigger impact. The right panel of fig. 
[12] shows that SM23 must decrease by about 1 keV in the 
INTEGRAL/DAMA-allowed region to compensate the 
effect of making €3 arbitrarily small. 



V. POSITRON PRODUCTION RATE AND 
ANGULAR PROFILE 

A. Rate from inelastic collisions 

The most recent determination of the observed 
positron annihilation rate in the bulge is 1.1 x 10 43 /s 
@. This value depends upon the assumed distance be- 
tween the sun and the galactic center; consistently with 
§ we take r Q = 8.5 kpc 41 1. For the predicted rate, we 
have 



Rp+ — — 



Ytot 



5 1 \ r 

x 



(28) 



where a is the cross section for XsXs — > XuXu (recall that 
Xs,u are the stable and unstable excited states). The 1/2 
is to avoid double-counting, and the abundance factor 
Ys/ Ytot is given by ([25]) . We integrate over a region of 
radius 1.5 kpc, corresponding to an angular diameter of 
approximately 10°. The observed profile, fig. M suggests 
that the signal falls below the sensitivity of INTEGRAL 
near this angle. 

The phase space average of av is given by 



(<rv) = / d 3 V!d 3 V2 f(vi) f(v 2 )cr\vi -#2! 



(29) 



We take a Maxwellian velocity distribution 

f(v) = 7Ve^>o 9 (v - v csc ) 
cut off at the escape velocity 

v 2 esc (r)=2v 2 (r) [2.39 + ln(10 kpc/r)] 
and having velocity dispersion 

Vo(r) 3 cx r x p(r) 



(30) 



(31) 



(32) 



with x = 1-64, and the normalization such that Vq{tq) is 
220-230 km/s. This form of vq is suggested by A-body 
simulations that include the effects of baryonic contrac- 
tion [29[ . Our choice of v csc follows ref. [42| ; see appendix 
C of that paper. 

In our previous work, the major challenge was to com- 
pute a since we were concerned with DM at the TeV 
scale, implying gauge couplings a g larger than the aver- 
age DM velocity v. In this case a nonperturbative cal- 
culation of a was necessary, since multiple gauge boson 
exchanges occur when v < a g , similarly to the Sommer- 
feld enhancement in DM annihilation (36[. However in 
the present situation a g <C v and a perturbative treat- 
ment suffices. 

We define some kinematic variables to facilitate the 
presentation of the cross section: 



5M> 



2:! 



(33) 



where SM23 is the small splitting between the two excited 
states, v is the DM velocity in the center of mass frame, 
and vt is the threshold velocity for X2X2 — > X3X3 excita- 
tions. The cross section for excitations can be expressed 
as 



v t Ml 

a t v rcl = 47ra;VA - 



where 



2rj 1 — 77 



(34) 



V = M 2 v 2 (2A- 1) 



9 

Mi; 



V = 2 



M^ t VA(A-l) 



x -rv— -i ■ ■, - v 

(35) 

Notice that cr+i; re i vanishes at threshold, A = 1. The 
related cross section for X3X3 X2X2 de-excitations is 



-v t M. 



where 



1-fj 



„ In — 

2 2fj 1 



V 



(36) 



2 2/ A ^ 2 MVVA(A + 1) 

P-M^ 2 (2A + l)+/i 2 ; fj = 2 x tV ^ 

(37) 

As expected, cr^v Te i can be obtained from cr^v Te i by chang- 

2 and 



ing 5M 2 3 — > —5M 2 3, which implies vf —> 
A — > — A (notice that UtyCTT = \/v 2 . . .). In the limit 
A — > 0, and substituting 8M23 ~ > 5M± S and /ii — > fii, we 
recover the zero- velocity cross section for downscattering 
through the large mass gap, 
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M (GeV) M (GeV) 

FIG. 12: To illustrate the effect of £2,3 on kinetic decoupling and the relic density of the excited state, left: contours of 
log R e + /Robs in the exothermic model, varying £2 between 10 _1 and 10~ 10 . ^12 is fixed at 250 MeV and other parameters 
are as in fig. [2] Right: similar plot for the endothermic model, with /i2 = 500 MeV. For clarity only the contours with 
log7? e + /R 00 s = are shown. Dependence on £3 is saturated for £3 > 10 -3 or £3 < 10 -5 . 



B. Rate from decaying DM 

We consider the scenario where the "unstable" state 
Xu is so long lived that it is already present in the galaxy 
due to its relic density, and decays with a lifetime t u 
greater than the age of the universe. Assuming that the 
511 keV gamma rays observed by INTEGRAL come from 
a central region of radius r c , the rate of positrons is 



R,+ = 



An 



M X T U JQ 



4ttC 



drr 2 p u (r) 
p kpc 3 ( Y,,, 



M x t u \Y t 



(38) 



where p u is the density of \ u and Y u /Y to t is the abun- 
dance of Xu relative to the total DM population. We 
assume the Einasto profile to obtain the dimensionless 
factor 
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FIG. 13: Contours of f, defined in (J39J), in the plane of the 
Einasto halo parameters, with r c /rQ — 0.176. £ is related to 
the volume integral of the DM density in the region of the 
INTEGRAL 511 keV signal, eq. gO}. 



( e (2/a)(r /r s )° 



\kpc / 
1 /ay 



r|2)-rfi.»(S 

a I \a a \r s 



(39) 



Matching R e + to the observed rate 3.4 x 10 42 /s, we find 
that the lifetime of Xu relative to the age of the universe 
(here we define tjj = 10 10 y rather than the actual age of 
the universe) 



2.1 x 10 5 



3Fn 
^tot 



5 GeV 



(40) 



To see how small £2,3 this corresponds to, we can 
rescale the bound (|5T|) to be derived below from demand- 
ing that similar decays of the "stable" excited state must 
take longer than t\j. It implies that 



c 2,3 



Yl\2 



S (7 x 10~") 



10" 



"0 



(t 



M2,3 



Ytot 
3KT 



5 GeV 



GeV J 

( P° 



0.1 \ vy 



SM 



Is- 



(41) 



where p© = 0.3 GeV/cm . The factor £ is plotted over 
a wide range of Einasto parameters in fig. 1131 showing 
that it is between 25 and 75 for reasonable profiles. With 
Y u /Y tot ~ 1/3 this gives t u ~ 10 5 times Ty. 



where 8M\ S - = 8Mi s — 2m e is the energy available for 
the decay. It is theoretically easy to achieve the desired 
rate of positron creation just by adjusting e 2 ,3 to this 
small value, since there is no other constraint on £2.3. 
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C. Angular distribution 



In this section we elaborate on the angular profile of 
the 511 keV signal in the case of scatterings only, since 
only there is it definitely necessary to consider the effects 
of positron propagation. The intensity of the signal as a 
function of angle is found by computing the line-of-sight 
integral (|10[) where the line is oriented along the direction 
x specified by angles 9, <j> relative to the galactic center. 
This expression assumes that the positrons decay at the 
same position where they were created. To model the 
effects of propagation before decay, we smear the angular 
distribution by averaging x weighted by some function 

/(COS0), 



I e +(x') = / dQf{x-x')I e +(x) 



(42) 



The integral over solid angle can be combined with the 
integral over the line of sight and rewritten in terms 
of a volume integral, dxdfl = d 3 x/x 2 , with the ori- 
gin of coordinates at the earth. Now, given that p 2 is 
strongly peaked near the galactic center, we can write 
p 2 ~ 5^\x — Xq), where xo is the position of the galactic 
center. Then we find that 



J e+ (0)~/(cos0) 



(43) 



The intensity has the same shape as the smoothing func- 
tion. As argued in section III Dt this is a good approxi- 
mation for the DM halo profiles that we are considering 
for the inelastic scattering mechanism. 

It is interesting to notice that even though p 2 looks 
like a delta function with respect to the measure d 3 x/x 2 , 
not so for the usual volume measure d 3 x. Indeed, the 
function r 2 p 2 (r) has a maximum near r even m 

the limit a —> 0. Therefore the total rate of positron 
production in the galaxy gets significant contributions 
away from the galactic center, although these are not 
counted in the observations of the bulge component upon 
which we focus in this paper, since only near the center 
is the intensity high enough to be detected. 



D. Regions of positron annihilation 

In the above discussion we have assumed that positrons 
are able to migrate to the regions where positronium 
forms and where they can subsequently annihilate. It 
is known from fitting the observed 7 ray spectrum that 
~ 97% of the positrons indeed form positronium be- 
fore annihilating [43l . |44| . This is because orthopositron- 
ium decays to three photons, and comparison of the 511 
keV line flux with the continuum level is consistent with 
nearly all annihilations coming from positronium rather 
than positrons encountering free electrons. The spectral 
shape also shows that most annihilations take place in 
warm (~ 8000 K) [45| regions, which may be mostly ion- 
ized [43| or else a combination of neutral and ionized 
regions [44] . 
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FIG. 14: Superposition of molecular gas regions of [Hj (green 
ellipses in center) on the intensity map of the INTEGRAL 
511 keV observations from ref. [g. Innermost ellipse is the 
CMZ (central molecular zone); outer tilted ellipse is the holed 
galactic bulge disk. 



Efforts have been made to independently map out the 
positions of the warm regions in the galactic bulge (GB); 
doing so could provide a consistency check on the above 
determinations, since then the morphology of the INTE- 
GRAL detection of the galactic bulge 511 keV gamma 
ray s should match the position of the warm regions. Ref. 
[46j | has modeled the spatial distribution of molecular gas 
in the GB based on CO emission data 47] for the central 
molecular region (CMZ) in the inner 150 pc, and borrow- 
ing an older model [48j based on Hi observations for the 
"holed GB disk" region extending to radii of ~ 1 kpc. 
It has been suggested that the warm neutral or ionized 
regions relevant for positronium ionization coincide with 
these molecular gas clouds [J]. 8 In figure [14] we have 
tranposed an outline of the CMZ and holed disk regions 
(fig. 4 of ref. HI) on the most recent INTEGRAL 511 
keV intensity map @. 

From fig. M it is clear that the INTEGRAL/SPI in- 
strument does not have sufficient spatial resolution to 
test whether positron annihilation really comes from the 
molecular gas regions. Furthermore, the assumption that 
these regions coincide with the warm ionized or neutral 
+ ionized regions of positronium annihilation is question- 
able. The direct measurements of the ionized component 
are based upon pulsar observations [50], which suggest 
the existence of warm H + regions with similar morphol- 
ogy to the molecular gas. But this is not considered to 
be a very reliable measurement of the ionized gas den- 



8 Ref. [5| assumes that positrons from the radioactive decays of 
supernova ejecta can be transported from the galactic disk into 
the GB to account for the observed 511 keV signal; however the 
validity of their model of electron transport has been questioned 
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sity in the GB due to the scarcity of pulsars in this region 
[5l| . Thus we do not know with a high level of confidence 
where the warm regions of positronium annihilation are 
really located. It is possible that they extend beyond the 
molecular gas regions identified by [461 ] . 

Because of the lack of very reliable information as to 
the spatial distribution of the warm ionized (+ neutral) 
regions, an uncertainty that is acknowledged in careful 
studies such as [46| , it is possible that decays of DM lead 
to positronium production in the vicinity of the initially 
produced positrons, so that the INTEGRAL signal could 
be a reflection of the underlying DM distribution. On the 
other hand if the positrons result from DM scattering, 
we have shown that they are initially produced within 
1° ~ 150 pc of the galactic center (dotted red curve of 
fig- El), and then transport of the positrons to larger radii 
is probably necessary to be consistent with the observed 
extent of the 511 keV signal. Positron transport in the 
galactic center has been extensively studied, and shown 
to depend sensitively on the largely unconstrained nature 
of the magnetic fields in this region. 

Ref. recently showed that, under reasonable as- 
sumptions, positrons can travel well outside of the GC 
before annihilating. If this is the case, the initially highly 
localized source from DM scattering will be widened to 
fill the interstellar medium. More information will be 
needed to attach a firm interpretation to the angular dis- 
tributions of the DM decays or annihilations. Obser- 
vations of the 511 keV gamma rays using a future in- 
strument with better spatial resolution would clearly be 
desirable for helping to settle these questions. In partic- 
ular, if a new measurement revealed stronger localization 
of the GB component of the 511 keV emission toward the 
galactic center, it would favor the DM explanation over 
astrophysical sources. 



VI. DIRECT DETECTION RATES 

In our computation of the 511 keV rate, we fixed the 
value of gauge kinetic mixing parameter e\ so as to match 
the direct detection rates determined respectively by ref- 
erences [IH and 24 1 for the exothermic and endothermic 
DM models. Although ei does not directly affect the rate 
of X2X2 X3X3 transitions, it does so indirectly, because 
of its influence (through kinetic decoupling) on the relic 
density of the stable excited state. Here we give details 
on the determination of e± in these two cases. We note 
briefly that the nuclear form factor is trivial for collisions 
studied here, so we will ignore it. 



DAMA transitions: 



OVi,cl 



4ttA 4 



(44) 



where p n — m n M x / (m n + M x ) is the reduced mass and 
A = 340 GeV. 9 In our model, the coupling is to protons 
only, and the analogous quantity is given by 



1 a 2 Mm 

a p>e i = lb7re 1 aa ff — 
Mi 



(45) 



To determine the value of t\ needed to match the ob- 
served rates, we must account for the coupling to protons 
only since the rate is proportional to A(Z f p + (A— Z) f n ) 2 
for atomic number and mass Z, A and relative strengths 
of couplings to protons and neutrons / p , /„ . In ref. [18| , 
the couplings were assumed to be f v = f n = 1, but we 
have f p = 1, f n = 0. Moreover, we have a different local 
density of the excited state than that assumed by [l8j 
because of the abundance factor Y s /3Yi (which also ap- 
pears in the positron rate ([2"51) ). and because we allow 
the local DM density to vary with respect to the fiducial 
value po — 0.3 GeV/cm 3 . The result is 



ei 



8ttA 2 Z 



1 Y tot po 
aa g Y s po 



1/2 



(46) 



If £2,3 <C ei so that ei determines the kinetic decoupling 
temperature of the DM, then Y s depends implicitly on 
t\ and (|46p must be solved numerically. The factor with 
A/Z depends upon which nucleus we are talking about, 
and is given by 2.4 and 2.28 respectively for I and Na. 
As [iH notes, scattering from Na nuclei is preferentially 
detected in our region of parameter space, so we choose 
the latter number. 

In figure [T5] (left panel) we plot contours of 6\ corre- 
sponding to the P2 = 500 MeV example shown in fig.[U to 
give a sense for how large t\ must be. Near M x — 4 GeV, 
5M 23 = 4.5 keV, where the INTEGRAL and DAMA 
rates best fit simultaneously, t\ ~ 10~ 327 , significa ntly 
larger than the generic estimate 10~ 5 given in ref. (18f . 
This is due to the A/Z correction, the fact that we need 
pi to be heavier than the nominal 100 MeV value as- 
sumed in [l8j], and that Y s /Y tot can be significantly less 
than 1 in our model. 



B. Endothermic dark matter 

Ref. :24]] finds that DM with a mass of approximately 
7 GeV and cross section on nucleons 



Exothermic dark matter 



cr„ = 2 x 10^ 4 pb 



(47) 



Ref. [18j determined the elastic limit of the DM- 
nucleon cross section needed to get the right rate of 



9 Note that taking this elastic limit negates the need to average 
over DM speeds as in 1 1491 . 
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FIG. 15: Contours of log ei (dashed) and (solid, in MeV) for [12 = 500 MeV, in the exothermic (left) and endothermic 
(right) models. Thick curve labeled "511 keV rate" is the contour where the predicted positron rate matches the INTEGRAL 
observation. Shaded areas are the DAMA-allowed regions. 



can simultaneously fit the DAM A/LIBRA annual modu- 
lation and the CoGeNT low-recoil events. Their allowed 
regions of a n versus M x are reproduced in figure [5] The 
logic for matching our cross section to theirs is similar to 
the exothermic case, except for the fact that endothermic 
scatterings are kinematically blocked if the DM velocity 
is below the threshold value 



v t = y/2SM 23 /nN (48) 

where /ijv = m^M x j '(mjv + M x ) is the nucleus-DM re- 
duced mass. We take this into account by doing the phase 
space average of av. The phase space factor in av that is 
sensitive to the threshold is \/v 2 — v 2 . For elastic scat- 
tering, this factor would be v. Therefore we match the 
quantity (|45|) that also appears in our slightly inelastic 
cross section to (jT7|) using 

<"> ( A ) 2po3Yl -... (49) 
where the averages over velocity are performed with the 

2/2 

Maxwellian distribution function / = Ne~ v /v o cut off at 
the escape velocity v osc . Since we are comparing with ref. 
|2 ij . we use their values vq = 230 km/s and v csc = 600 
km/s for this part. Once again, scatterings from Na are 
preferentially detected, so A/Z = 2.28. Similarly, we 
take the threshold velocity for sodium in the above. 

Figure [T5] (right panel) shows contours of ei for the case 
of ^ 2 = 500 MeV. In the overlap region for INTEGRAL 
and DAM A, e\ = 10~ 5 . This is smaller than required in 
the exothermic model because the corresponding value 
of fj,i is smaller, and also the cross section ([T7]) is ap- 
proximately 0.15 that in (l44l) . Kinetic equilibrium of the 
DM with the SM in the early universe is not as efficiently 
maintained by \2 ^ X3 transitions in this case. This is 
why the relic density of \i is sensitive to the value of 



the other nonvanishing kinetic mixing parameter for the 
endothermic model, whereas it is practically insensitive 
in the exothermic case. 



VII. ASTROPHYSICAL CONSTRAINTS 

In this section we address the astrophysical and cos- 
mological constraints on our proposal that are comple- 
mentary to the 511 keV and direct detection signals, as 
well as to laboratory constraints from electron beams. 

A. Lifetime of metastable state 

We need to insure that the "stable" excited state \s 
is either truly stable or else sufficiently long-lived. The 
most dangerous process is Xs ~^ Xi e+e - At the phe- 
nomenological level, we suppress this by setting the ki- 
netic mixing of the gauge boson that mediates this pro- 
cess to a sufficiently small value. The rate for this de- 
cay can be computed analytically with the approxima- 
tion that 5 Mis- — SMis — 2m e is sufficiently small for 
the final state particles to be nonrelativistic. Then 

T Xa = 2a g ae 2 m 2 e 5M? S _ ^ 4 (50) 
Demanding that r s exceed 10 10 y requires that 

\ a g J VI GeV/ V, <5M ls _ J 

(51) 

In the endothermic model, we can set 63 = at tree 
level by removing the A 3 Higgs boson, but there seems 
to be no symmetry to ensure that €3 is not generated 
by loops if the other two mixing parameters are nonzero. 
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Nonetheless, we are not able to find an example of a loop 
diagram that generates nonzero any that superficially 
seem promising vanish because of Furry's theorem. In- 
stead, we find a one-loop process where B% acquires a 
magnetic moment coupling to the electron, fisB^ea^e. 
The decay of \2 proceeds by B3 exchange in the one-loop 
diagrams of fig. 1161 These diagrams would cancel exactly 
if /i2 = Hi, so the magnetic moment can be estimated as 



(/'2 



47T^ 2 



^) ln A 



(52) 



where p, = \{^\ + ^ 2 ) and A denotes the hidden SU(2) 
symmetry breaking scale, above which the kinetic mix- 
ing of -B1.2 is replaced by an interaction with the triplet 
Higgs fields. The squared matrix element of fig. [16] can 
be estimated as \M\ 2 — g 2 fi 2 B Al 2 mlSM 2 2 fi^ 4 , and the 
decay rate in the limit of small SM12- = SMi 2 — 2m e is 



r 



X2->Xie"'"e- 



^-li B m\ 8M 2 2 SM 2 2 _^ 



(53) 



For a g ~ 10 4 and e\ ~ e 2 ~ 10 3 , <5/i ~ /13 ~ fj, ~ 1 
GeV, SM 12 - ~ 0.1 MeV, A - 10 GeV, we find a lifetime 
of 10 26 s, much larger than the age of the universe. There- 
fore it seems technically natural to neglect the dangerous 
kinetic mixing term and assume the "stable" state is suf- 
ficiently long-lived. As it turns out, a careful calculation 
is even more suppressed; see appendix [Q 





FIG. 16: Decay of metastable X2 state due to nonzero £1,2- 

The X3 ~^ X2X decay in the exothermic model is not 
problematic, since the mass splitting is much smaller and 
the only available decay channels are with X = 7, X — 
vVj. an d X — 37. These have been studied previously 
[lIlHHHH- The single-photon decay has a lifetime longer 
than the age of the universe for the value of e± required 
for direct detection; since it could be observed, we discuss 
it in more detail in section IVII Bl below. 

The partial width for the vv final state is easy to es- 
timate in analogy with (f50|) . For this channel, there is 
an additional suppession in the kinetic mixing. First, 
the mixing of B\ with the Z boson current has an ex- 
tra factor of n\/m 2 z [Hj], and, second, the SM Z boson 
mixes with the B\ current with opposite sign such that 
the two mixings nearly cancel at small energy-momentum 
transfer (see appendix [B] for discussion of the same can- 
cellation in \ v scattering). The ensuing bound on e± is 
much weaker than that on e 2 ; practically speaking there 
is no constraint. 



The decay X3 — > Xi + 37 is due to the operator 
~ (ea 2 /Q0m 4 )B^ v ' F^F 2 induced by an electron loop, 
similar to the Euler-Heisenberg F 4 interaction in QED. 
(Furry's theorem forbids a term of the form B\F 2 and 
B\F mixing is already taken into account by diagonaliz- 
ing the kinetic terms.) The rate is suppressed by SM 2 ^ 
[l3L I22I ] , leading to lifetimes that far exceed the age of 
the universe for the small ~ 5 keV splittings relevant to 
our exothermic model. 



B. Single-photon decays of metastable state 

The exothermic version of our proposal faces the chal- 
lenge that the excited state can decay by emission of a 
single photon, via X3 ~ ^ r X2l- The origin of this de- 
cay was pointed out in [35|]: the nonabelian term in 
the field strength B^ v leads to interactions of the form 
eiB^B^F^ with the photon, from the gauge kinetic mix- 
ing operators ([1]). This can be put into a loop diagram 
which results in a transition magnetic moment X2-X3 1 



M23 = 128tt 2 M x 



l 



(54) 



where /i is of order ^2,^3 • Therefore, there is a decay 
channel X3 ~^ X27- The rate is 



#^M 23 ) 2 



(55) 



For ei 10~ 5 , M x 4.5 GeV and 5M 23 = 5 keV, the 
lifetime is 4 x 10 19 s which is much longer than the age 
of the universe. However, this is not necessarily enough 
because such photons could be visible in astronomical 
searches. 

Of the various instruments that could be sensitive to 
low-energy 7 rays, INTEGRAL/SPI comes the closest. 
Ref. [HH gives limits on the intensity of gamma ray lines 
that could come from such decays in the galaxy; however 
INTEGRAL'S sensitivity cuts out below 20 keV, making 
our scenario just out of reach. Interestingly limits on 
the diffuse gamma ray background put a (5M23-dependent 
lower limit on the partial lifetime times the mass [53] of 
approximately 



r 7 M x > 1 X UT 



3*3 

Ytat 



8M- 



2:i 



10 keV 



1.2 



GeV s (56) 



for SM 23 > 10 keV. Again because of INTEGRAL'S en- 
ergy sensitivity, data is not given for lower photon ener- 
gies. Nevertheless, extrapolating the bound to SM 23 = 5 
keV gives 4.3 x 10 19 GeV s, which is not even five times 
less than our nominal value 1.8 x 10 20 GeV s, assum- 
ing M x = 4.5 GeV. Therefore an instrument sensitive to 
these lower energies might detect this low-energy photon, 
which is in the x-ray part of the spectrum. 

In fact, observations of the galactic center by the Chan- 
dra x-ray telescope [54| may rule out this particular 
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the same form, except for the substitutions of e and /i 
by the corresponding quantities for B3, in the exother- 
mic model; for the endothermic model, ([50)) applies as 
written, to the unstable excited state. Laboratory ex- 
periments constrain the appropriate e$ to be ;$ 10~ 3 , so 
the lifetime could be 10 10 times shorter then the above 
estimate using e ~ 10~ 8 , thus on order of 1 y. This as- 
sumes the large mass splitting is only 1.1 MeV. With a 
2.1 MeV splitting one gains a factor of 10 5 in the rate 
due to the larger phase space, giving a lifetime of several 
hundred seconds. It cannot be significantly smaller in 
our model. 



FIG. 17: Chandra spectrum from region near galactic center, 
where ~ 5 keV x-ray from \3 ~* X27 decay might manifest 
itself. 



model. Observations are presented for a region of size 
35 arcmin 2 that is 7.5 arcmin away from the GC. No ev- 
idence of an unidentified line is observed in the 1 — 8 keV 
band (fig. I17[) , and the continuum seen there is modeled 
by thermal sources with a flux of6x 10~ 12 erg cm~ 2 s _1 . 
We can compute the expected flux by integrating over the 
line of sight and the solid angle (dSl — dipdcosip) of the 
observed region [53j], 



D. Single-photon decays of unstable state 
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where I = e^/^y" J d cos ip J dl e- 1 - 2 / ^ +l -Uv^^Y , 

y = rQ/r s and I = l/r s . We numerically integrate over 
an annular region of similar area and displacement from 
the GC to the observed one, using the Einasto param- 
eters ((7|) to find I = 10~ 4 . Using the value of M x t 1 
determined above, this gives F t h = 0.001 photons cm~ 2 
s _1 (for the maximal ratio ia/ltot' The corresponding 
energy flux for a 5 keV mass difference is 10~ n erg 



The decay mechanism discussed in section IVIIBI was 
originally conceived for the decay of the unstable state in 
ref. [Hj]. This goes through the MeV-scale mass gap, so 
the photon in this case is a gamma ray. For definiteness 
let us consider the endothermic model, so X3 is the un- 
stable state and the relevant decay is X3 — > xil y i a the 
/U13 transition magnetic moment, which is proportional 
to £2 7 in analogy to The partial decay rate is the 

obvious generalization of (1551) . The branching ratio for 



the single photon decay relative to that into e + e is [35 1 



ai/a 



7 8192tt 2 M?(dM 13 _) 3 (SM 13+ f 



In 2 (58) 
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where /1 ~ Mi, ^3) $M\s± = <5A/ 13 ± 2m e and e = 
2.71828 .... The resulting photon might be observed by 
INTEGRAL in the diffuse 7 ray background. In ref. (35j . 
a bound was derived, which however overestimated the 
sensitivity of INTEGRAL to the signal. We therefore 
reconsider it here. 

The analysis of ref. [52[ is particularly relevant for us, 
since they searched for line sources from the galactic cen- 
ter region, having a spatial distribution similar to that of 
the 511 keV line. They limit the flux of such a line, for en- 



cur 2 s -1 , not quite 2 times greater than the observed ergies between 1 and 2 MeV, to less than -3x10 5 cm 



continuum flux. 

In the above estimate we did not take account of ab- 
sorption of the decay signal, which would help to soften 
the discrepancy, especially if the photon energy is some- 
what lower. It may also be possible to evade the problem 
by extending the gauge group to SU(2) x U(l) [55J and re- 
place the kinetic mixing of B\ by that of the extra U(l) 
gauge boson; this would remove the ^.23 transition mag- 
netic moment. Notice that this problem does not affect 
the endothermic model because X3 ~^ Xi e+e ~ proceeds 
much more quickly, as we discuss in section [VII CI below. 



This is to be compared to the flux from positrons, 
— 3.6 x 10~ 3 cm~ 2 s _1 . Therefore BR 7 should not be 
greater than about 0.01. However, using the typical val- 
ues of parameters of interest for our present application, 



we find BR 7 
searches. 



10 



-7 



far below the sensitivity of current 



E. 



Cosmic ray and CMB constraints 



Lifetime of unstable state 



In passing, we can also estimate the decay rate for 
the unstable excited state into e + e~ using (J5UJ. It has 



DM annihilations can occur even after freezeout, with 
the production of gamma rays or charged particles that 
can have an observable effect. In the class of models we 
consider, the DM annihilates directly into hidden sector 
gauge bosons, which in turn decay into any charged SM 
particles that are sufficiently light. Gamma rays emerge 
only as secondary products of these charged particles. 
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FIG. 18: Adapted from ref. [59j, showing constraints on 
the DM annihilation cross section versus mass from opti- 
cal depth ("CMB") and excess heating of the intergalatic 
medium ("gas"). Regions above the diagonal lines are ex- 
cluded. "Relic density" region indicates the desired value of 
the cross section for the correct thermal abundance. 

Their contribution to the diffuse gamma ray background 
can potentially give interesting constraints [5(| , but cur- 
rently the uncertainties from details of structure forma- 
tion do not allow one to draw firm conclusions. The 
production of antiprotons in the galaxy gives more def- 
inite constraints, which can be quite stringent (57j . To 
avoid them, we need to assume that the gauge bosons 
which mix with the photon are lighter than 2m p so that 
pp pairs are not produced. 

Charged particles that are produced around the time 
of recombination reionize the plasma and change the op- 
tical depth to the surface of last scattering, a quantity 
that affects the Doppler peaks of the cosmic microwave 
background [H, [59j] . The effect is particularly strong 
for DM with mass M x <, 10 GeV, as in the exothermic 
proposal for DAM A. Ref. [II] shows that such DM is 
marginally ruled out if it decays exclusively into e + e _ , 
while it is marginally allowed if it decays into heavier 
charged particles (which decay into electrons that are 
less energetic than if they were primary products) . The 
relevant bounds are reproduced in fig. [T5] In most of 
our examples, the gauge bosons can decay into muons 
and charged pions, so the branching ratio into electrons 
will be suppressed and the bound should be somewhere 
between the two cases shown in fig. [T5] Ref. [58[ also 
derives a bound from excess heating of the interstellar 
gas, which is more stringent than the CMB bound in 
this small M x region, but which is also more subject to 
uncertainty because of its dependence upon assumptions 
about the details of structure formation. 



F. Long-lived gauge bosons and nucleosynthesis 

It is interesting to consider possible effects of the hid- 
den sector gauge bosons in the early universe. Decays 
around the time of big bang nucleosythesis or later can 
be deleterious, although it is also possible to improve the 



predictions of big bang nucleosynthesis (BBN), notably 
for lithium [60']. In our models, Bi is usually the lightest 
gauge boson, and it couples to electrons (and muons and 
pions) with e\ ~ 10~ 5 , fixed by the rate of direct DM de- 
tection. Its decay rate is therefore of order aef/ii which 
for (j, i ~ 500 MeV leads to a lifetime of 10 -10 s, which is 
clearly harmless. 

On the other hand, if £2 ~ 10 -10 in the exothermic 
model, where we have the constraint ;$ 10~ 8 , then B2 
can have a lifetime greater than 1 s and possibly be rele- 
vant for nucleosynthesis. The question is whether its relic 
density is large enough to have an effect. We have com- 
puted the cross section for B2B2 — > B\B\ using FeynCalc 
jolj . The cross section as v — > can be expressed as 

7ra 2 

™ = Trifix) (59) 

where x = 1 — ^2/^3 = 2(5Af23/a 9 M3 and / has a mini- 
mum value of 18 at x = (treating (g]) as an equality to 
eliminate /ii). This neglects dark Higgs exchange in the 
s-channel, but we have checked that including it makes 
no dramatic change unless the virtual Higgs goes on shell. 
For typical values we find that the standard relic abun- 
dance calculation gives a freezeout temperature around 8 
MeV for B2, and an abundance 10 -2 times smaller than 
that of the baryon asymmetry. This is too small to have 
any effect on BBN. 

Ref. [gFj] point out that a more likely candidate for 
giving interesting effects is the dark Higgs bosons. In 
particular, if there exist a Higgs boson that is lighter 
than the gauge bosons, it would decay into 4 leptons 
through emission of two virtual gauge bosons, with a rate 
suppressed by e 2 a g a 2 (m/ l //i) 8 . The annihilation cross 
section is suppressed for similar reasons. This can more 
naturally give long-lived relics (on the time scale of BBN) 
that could solve the lithium problem. 



VIII. LABORATORY SEARCHES 
A. Beam dump experiments 

An interesting feature of the class of models we con- 
sider is that they can be tested in proposed low-energy 
laboratory experiments. A beam dump on an absorbing 
target can produce the weakly interacting B bosons that 
mix with the photon. These can decay into e + e~ or other 
charged particles before reaching the detector, providing 
a signal not present in the standard model. 

In our scenario, two of the three colors of bosons should 
mix with the photon: Bi, with strength e\ ~ 10~ 3 — 10~ 5 
to get the right rate of direct detection, and either B 3 
or B 2 , depending upon whether \i or X3 is the stable 
excited state. Let us denote the corresponding mixing 
parameter by £3,2- We noted in section IIV CI (figure 
1121) that this parameter is essentially unconstrained. If 
£3,2 ~ 10~ 6 then the effects of B3 2 could be discovered 
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in laboratory searches. But since we have more definite 
predictions for B\, we will focus here on its discovery 
potential. Moreover we have argued that there are cer- 
tain advantages to having very small values of 63,2 which 
could make laboratory detection of B3 2 impossible for 
the present. 

The authors of ref. 62] has recently summarized 
the current experimental constraints in the plane 
(where fa is the mass of the relevant gauge boson), and 
they have also proposed strategies for new experiments 
that can cover more of the still-allowed region in this 
plane. Fig. 1191 reproduces some of their results. On top 
of these we plot several examples of predictions from our 
cndothermic and exothermic models (circles containing 
"n" or "x" respectively), corresponding to those shown 
in figures [5] (see also fig. [To]) . Almost all of these 
points are contained within the contours denoting the 
reach of feasible new experiments suggested by ref. [62j : 
the solid (blue) line denoting the high resolution, high 
rate trident spectromenter, and the dashed (red) one for 
the thin-target with double arm spectrometer. It is sug- 
gested that such experiments would be feasible at several 
existing laboratories, including Jlab (Thomas Jefferson 
National Accelerator Facility), SLAC (Stanford Linear 
Accelerator Center), ELSA (Electron Stretcher and Ac- 
celerator), and MAMI (Mainzer Mikrotron) . We see that 
only one of our examples (the right-most "x") would lie 
outside of the reach of the proposed experiments. This 
corresponds to the extreme case where /i2 = 2 GeV in 
fig. [5] The more typical models would therefore be in 
the discoverable region. 



B. Invisible width of Z boson 

The nonabelian gauge kinetic mixing portal ([JJ pro- 
vides two invisible decay channels for the Z boson: Z — > 
Aj.B, since e, stands for the VEV of the Higgs triplet 
A, over the heavy scale A (eq. (fTTj) '). and Z — > BjB^ 
where k are a cyclic permutation of 1,2,3. The lat- 
ter arises because the nonabelian field strength B^ v con- 
tains geijkBjB%. Considering the first process, the par- 
tial width is 



96ttA 2 



(60) 



in the approximation mz is much greater than the masses 
of the decay products. Demanding that this be less than 
the experimental error on the invisible Z width, 1.5 MeV 
63], we find that A > 1.3 TeV. For (A) - 10 GeV, this 
leads to the bound e ^ 10 -2 , which is less stringent than 
other constraints shown in fig. [19] For the Z — » BjB^ 
channel, the partial width is of order a g e 2 m z . This leads 
to a weaker bound on e than does the Z — > AiBi channel. 
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FIG. 19: Potential for discovery of light mixed gauge bosons in 
plane of kinetic mixing parameter e and gauge boson mass ^t. 
Shaded regions are ruled out by existing laboratory or astro- 
physical constraints. Unshaded enclosed regions denote the 
reach of experimental strategies proposed in ref. (621 ] . Circles 
containing "n" or "x" are typical predictions of our endother- 
mic or exothermic DM models, respectively. 



IX. CONCLUSIONS 

If the anomalous 511 keV gamma rays from the galac- 
tic center are truly distributed in an axisymmetric man- 
ner, as suggested by the INTEGRAL observations, this 
provides strong motivation to seriously consider DM de- 
cays or scatterings as their source, rather than localized 
sources such as supernovae or x-ray binaries. A new mea- 
surement with higher spatial resolution would be very 
desirable to help settle this question. In the meantime, 
it seems worthwhile to explore possible DM interpreta- 
tions, especially if they can explain more than just the 
511 keV signal. In the present work we have shown how 
a three-component DM model with two mass splittings 
and a hidden SU(2) gauge boson might address both the 
511 keV observation and indications of DM detection by 
DAMA/LIBRA and possibly CoGeNT. 

The scenarios we have presented involve slightly inelas- 
tic nuclear scatterings, either endothermic or exothermic, 
in the direct detection experiments: X2,3^ — > X3,2^'- 
The endothermic version with M x = 4 GeV is under 
stronger pressure from astrophysical constraints from the 
CMB ffigure lTS]) and especially from the decay X3 - > X27, 
not observed by Chandra (section I VII Bp . The latter 
could possibly be softened by some modification of the 
particle physics model, such as extending the gauge group 
to SU(2)xU(l). The exothermic model also requires a 
more cuspy halo than the endothermic one to get the 
observed 511 keV rate from X3X3 ~^ X2X2 scattering, al- 
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though still consistent with examples from TV-body sim- 
ulations that take into account compression by baryons. 

We have highlighted two distinct mechanisms for get- 
ting the 511 keV signal: either inelastic X2,3X2,3 - ► 
X3,2X3,2 scatterings followed by X3,2 — > Xi e+e ~ decays, 
or the decay process by itself when X3,2 has a lifetime 
of order 10 5 times the age of the universe. Whereas the 
first mechanism requires some mutual adjustments of the 
particle physics and DM halo parameters to get the right 
rate, the second is more easily arranged by just tuning 
the gauge kinetic mixing parameter £2,3 ~ 10 -11 that 
controls the decay rate. The decay mechanism points 
to the exciting possibility that the angular profile of the 
511 keV signal is actually a picture of the DM halo pro- 
file in the inner part of the galaxy, if positron diffusion 
is a negligible effect. The scattering mechanism on the 
other hand requires significant positron diffusion, or else 
propagation of the excited DM state before decay, since 
otherwise it predicts too narrow angular profile. It is in- 
teresting that our model can naturally explain such long- 
distance travel of the excited DM prior to its decay, by 
tuning £2,3 ~ 10 -7 . 

A very encouraging aspect of these proposals is their 
testability in low-energy electron beam dump experi- 
ments. The kinetic mixing parameter ei ~ 10~ 4 and the 
mass \x\ ;$ 1 GeV of the gauge boson mediating the di- 
rect detection scatterings are in prime territory for their 
discovery by such experiments, which could be done at 
existing laboratories. The models presented here are also 
potentially rich in consequences for cosmic rays, the dif- 
fuse x-ray or gamma-ray backgrounds, the CMB, and big 
bang nucleosynthesis. 
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out, so the gauge and Higgs bosons may be treated as 
massless. The effects of symmetry breaking are discussed 
at the end. For reference, we will consider a general gauge 
group and general representations for both DM and the 
Higgs. 

We consider first DM annihilation to Higgs bosons. 
This process is mediated by t- and u-channel diagrams 
involving Yukawa couplings at each vertex, s-channel di- 
agrams with an intermediate Higgs particle connecting 
a Yukawa coupling at one vertex to a scalar potential 
vertex, and an s-channel diagram with an intermediate 
gauge boson coupling to the DM and Higgs particles at 
either end. Since we are concerned in this paper with ei- 
ther parametrically small or vanishing Yukawa couplings, 
we assume that the s-channel diagram with an interme- 
diate gauge boson dominates. 

Consider incoming DM states Xi m representation R 
and outgoing Higgs states A/ in representation R'; the 
incoming momenta are pi and outgoing momenta q\. The 
matrix element is 10 
,2 



M = 



if! 



(Al) 



(Pi+Pj) 2 

Once summed over outgoing colors and averaged over 
incoming colors and spins, it is 

|2 



1 9" 



■Id 



tr 



R ■ 



MxKtt-iMi-MxKtj-ij) 



rjib rjiCl rjiQ rjiCl 

X 1 ji 1 ij 1 JI 1 IJ 



(A2) 



where dn is the dimension of representation R. The color 
sums both take the form 

ti R T b T a = {d R /d adj )C 2 {R)5 ab , (A3) 

where C 2 is the quadratic Casimir; the sum over the ad- 
joint indices gives S ah 8 ab — d a dj . For nonrelativistic dark 
matter at center-of-mass velocity v — v rc \/2 and scatter- 
ing angle 6, the amplitude becomes 



\M\ 2 = 



dadjd 



C 2 (R)C 2 (R') 



±{1 
2 V 



v 2 cos 2 ( 



(A4) 

Annihilation to gauge bosons (of colors a, b and mo- 
menta q a ,b) receives contributions from s-, t-, and u- 
channels. The amplitudes for each channel are (as in 
351) 



Appendix A: Annihilation amplitudes and rates 

In this appendix, we derive the invariant amplitudes 
(squared) for annihilation of DM particles x to both 
gauge and Higgs bosons, including the lowest order cor- 
rections due to dark matter velocity, which we use to 
find the relation between the thermal relic density and 
the dark gauge coupling a g in section II V Al Including 
annihilation to Higgs bosons extends and corrects the re- 
sults listed in (35[; in addition, we correct the final state 
polarization and color sums carried out in (35|. To keep 
the final result simple, we will first assume that symmetry 
breaking occurs at a lower temperature than DM freeze- 



= ?-v j T? i j x u i r bc £l(a)EUb) [r,^(q b - q a ) x 



M t =1 



M u =» 



Ml 



Ml 



10 We use a mostly plus metric with Dirac algebra {7^,7"} 
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We need to account for both direct squares and cross 
terms in the amplitudes. After some tedious algebra in- 
cluding Dirac traces, we find the following results for the 
(color and spin summed and averaged) square amplitude: 

\M S \ 2 = ^tv R T c T d f abc f abd (~ + ^v 2 (ll 



-5cos26>) , 



\M t 



2_l 

4 



tr fl (T a T b T b T a ) {l±vcosB + v 2 ) 



M s Mt = l -^-f abc tT R (T c T b T a ) (2 - v 2 si 

"'R 



2 sin 2 . 



M S M U = - l -§rf abc tT R (T c T a T b ) (2 - v 2 sin 2 ■ 

"-R 



2g 4 



MtM u = ^-tv R (T a T b T a T b )v 2 
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Here, we have expanded around zero DM velocity as for 
annihilation to Higgs bosons. 

We can evaluate the group theory factors using 
jabcpbd = c 2 (adj)S cd and T a T a = C 2 (R), antisymme- 
try of the structure constants, and the group algebra. We 
find 



t 



\M S \ 2 = ^C 2 {adj)C 2 {R)(~ + \v 2 [\\ 



-5 cos 20) 



\M t , u \ 2 = ^-C 2 {R) 2 (l ±v cos 6 + v 2 ) , 
dR 

M s Mt = ^—C 2 (ad 3 )C 2 (R) (2 - v 2 sin 2 9) , 

Mr 
4 

M S M U = ^—C 2 (adj)C 2 (R) (2 - v 2 sin 2 9) , 

MR 

M t Mu = 2jL (c 2 {R) 2 + ~C 2 (adj)C 2 (R)^ v 2 (A7) 



so we find 



\M\ 2 = C 2 (adj)C 2 (R) (I + v 2 (l- cos 2 6)) . (A9) 
m R 

As an example of this effect, we can consider X1X2 ~> 
BiB 2 scattering (with fixed colors for SU(2) triplet DM). 
Taking just the kinematical factors, the amplitude for 
annihilation into gauge bosons is given by \M S + M u \ 2 = 
5<? 4 /4, as in [35[. The ghosts subtract g 4 /8 for a total of 
\M\ 2 — 9<7 4 /8, in agreement with the massless limit of 
the amplitude in the symmetry breaking phase. 

We are primarily interested in SU(2) triplet DM with 
triplet and quintuplet Higgs fields, and we now specialize 
to that case, assuming N3 Higgs fields in the triplet and 
A^5 in the quintuplet. The total squared amplitude, with 
gauge and Higgs boson final states added incoherently, is 
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For identical nonrelativistic initial particles and identical 
massless final particles, the differential cross-section is 



da 
dtt 



1 1 



4a 



2 64tt 2 s \ fi 



\Mf 
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where the factor of 1/2 is due to overcounting identical 
final states in the color sum. (Here we have assumed 
that the Higgs particles are in a real representation; if 
not, drop the factor of 1/2 for the Higgs final states.) 

So far we have worked at tree-level and in the CM 
frame, but there are small corrections to both approxi- 
mations. First, the annihilation cross section experiences 
Sommerfeld enhancement since the DM is nonrelativis- 
tic. Under the assumption that the gauge symmetry is 
not yet broken, the boost factor is 



S(v) 



ira g /v 



1 



o — 7ra g /2v 
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In the above discussion, we took the sum over gauge 
boson polarization vectors to give the metric for simplic- 
ity; this includes longitudinal and timelike polarizations 
as well as the transverse ones. In a nonabelian gauge the- 
ory, the unphysical polarizations do not automatically 
vanish when contracted in the amplitudes, so we must 
correct for their inclusion. We can do this by subtracting 
the squared amplitude for ghost production. The ampli- 
tude is 



M - 9 -^ T M a u if abc 



(A8) 



In the parameter space appropriate to the thermal DM 
freezeout, Tra g /v < 1, so we treat it as a small parame- 
ter and expand S(v) ~ 1 + Tra g /2v. Next, since the CM 
frame is not quite the rest frame of the cosmic fluid, we 
must include the Lorentz transformation of the cross sec- 
tion. For nonrelativistic center-of-mass velocity V, this 
correction takes a — > a(l — V 2 sin 2 0/2), where <f> is the 
angle between V and v. 

Finally, we must average the cross section over the DM 
velocity distribution, which is Maxwell-Boltzmann. In 
terms of the V and v, this average takes the form 
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where the relative velocity v rc \ = 2v. We have ex- 
panded the result to first order in the small parameters 
cx g y/ 'M x /T and T/M x . In the sequel, we will treat the 
ratio M x /T as roughly independent of the DM mass and 
use an approximate value of M x /T ~ 20. 

We should also address the issue of symmetry breaking. 
We will assume that the gauge bosons are light compared 
to the DM and that Higgs bosons are either light com- 
pared to the DM or too heavy to be produced in DM 
annihilation. The annihilation amplitudes are slightly 
modified by the gauge boson masses (in a manner that 
does not respect the gauge symmetry, of course). For 
example, X1X2 B1B2 annihilation has 
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at zero velocity and lowest order in the gauge boson 
masses. This deviates from the leading order massless 
result only in as much as the gauge boson masses differ 
from each other. Note that the Goldstone boson states 
become longitudinal gauge bosons in gauges such as the 
unitary gauge. The significant effect of symmetry break- 
ing is to change the kinematical factors in the cross sec- 
tion. Using the average mass p. for all the light gauge 
bosons and Higgs states, 



|<f a | = M x (l + v 2 /2) J M*{1 + v*) - p? 



M v \ 1 



2 \ V2, 



M l 



M 2 v 2 
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For symmetry breaking at a small scale compared to the 
DM mass, this just renormalizes the cross section by a 
factor of (l-pf /M*) 1 / 2 . 

We can now compare the cross section to that required 
for the correct relic density of dark matter. Under the 
normal assumption that only SM particles are lighter 
than the dark matter, the required cross section is ap- 
proximately (av Ie \)o ~ 2.84 x 10~ 26 cm 3 /s, with a log- 
arithmic dependence on DM mass. We will take this 
central value. However, our dark matter models contain 
extra light degrees of freedom, which affects the required 
cross section in two ways. First, the Hubble parameter 
is greater at a given temperature, which causes freeze- 
out to occur earlier. Second, due to heating of photons 
by annihilation of these light degrees of freedom (see Big 
Bang nucleosynthesis constraints given in [35]), the uni- 
verse expands more between freezeout and the present 



day, which means that the freezeout density of DM must 
be higher than in minimal DM models. As a result, the 
desired cross section satisfies 



For DM masses near M x = 5 GeV, this is 
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v /l + (6 + 3A 3 + 57V 5 )/61.75 
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Including all these corrections, we can write the desired 
cross section in terms of the gauge coupling as 



(av lel ) = {Aa 2 + Ba 3 )/M 2 



which has the iterative solution 
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Our SU (2) models have 
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In the end, corrections due to the initial velocity of the 
DM particles contribute at the 5 — 15% level, while cor- 
rections from Sommerfeld enhancment contribute only 1 
part in 10 4 due to their additional dependence on a g 
(since the coefficient B is of the same order as A). 



Appendix B: Kinetic coupling to SM 

In this appendix, we find the freezeout temperature of 
the dark matter kinetic coupling to the Standard Model. 
For convenience calculating phase space factors, we con- 
sider downscattering xa,s/ — ^ Xifi labeling the Xj m0_ 
menta as pj and the initial and final energy and momenta 
of the SM fermions / as Eij and qij respectively (simi- 
larly for other subscripts). Here, the relevant initial DM 
state is \3 f° r t ne endothermic case and xi f° r exother- 
mic. In this appendix, we consider the endothermic case, 
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but it should be clear that these results apply equally 
well in both scenarios. At temperatures under 100 MeV 
(and, in particular near the important scale of <5Mi3 ~ 
MeV), only scattering from e^- will be important. Fur- 
thermore, if this process occurs roughly once per Hubble 
time for each of the more massive DM particles, it effi- 
ciently maintains the distribution of these two DM states 
given by the Boltzmann factor at the SM temperature. 
Also, if XX scattering is still in equilibrium, this reaction 
can insure that all the DM states maintain a thermal ve- 
locity distribution at the SM temperature. (As we have 
seen in the main text, XX scattering typically freezes out 



later than this process.) We note that the same calcu- 
lations apply for scattering between the two top states 
Xsf X2f with the appropriate replacements of £2, ^2, 
and 5M\3. As shown in figure [TTI the cross section is 
only slightly smaller for keV mass splittings. 

The matrix element for the scattering process shown 
in figure [10] is (taking electrons for specificity) 
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ffee 2 



(P3 - pi) 2 - 14 

After the spin sum and average 
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With the replacement m e — > m p , 5M 13 — > 5M 23 , and 
M2 - > fJ-i, this is also the result for inelastic DM scat- 
tering off protons, which is relevant to direct detection 
experiments. 

Since we are interested in temperatures much less than 
the DM mass M x (and all other energy scales are also 
much less than M x ), we can work in center-of- momentum 
(CM) frame up to overall errors of order X /T/M x <C 1 in 
the cross-section compared to the cosmic rest frame. To 
lowest order in T/M x , the final electron energy is 



Ef = E t + SM- 
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I q f | 2 =5M 2 3 + 2SM 13 E l + \q l \ 
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to lowest order. Due to cancellations in (|B2I) . we will 
need s to the same (second) order in small quantities: 



s = M 2 +2M X (SM 13 + S i )+5M 1 2 3 +2<5M 13 ^+m^+2|g i 



To lowest order in small quantities, we find 



[2 = 16 ( g ee 2 ) 2 M 2 
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Then the Mandelstam t ranges between t_ and t + satis- 
fying 



t± = -m\ 2 - 2SM 13 Ei 

±2\%\ [SM 2 3 + 25M 13 E t + \q t \ 2 } 1/2 
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for the (spin averaged and summed) squared amplitude. 
Since da/dt — (1/ 647rs) |7W | 2 /|<7i | 2 and the relative veloc- 
ity is dominated by the electron velocity, we find that 
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We are especially interested in whether the DM and 
SM can maintain kinetic equilibrium at temperatures 
T < m e , since those temperatures are relevant for x~ X 
downscattering. Then temperatures are parametrically 



less than /X2, so 
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times corrections of relative order t///| and a/|<&|/M x , 
where u re i is the CM frame electron speed. At a fixed 
temperature T, we find the total scattering rate by inte- 
grating over the Fermi-Dirac distribution for marginally 
rclativistic electrons (since T ~ m e ). Assuming a ther- 
mal origin for the dark matter abundance so a g oc M x , 
this total rate can be written in terms of a normalized 
rate T as 



in equation (IB7|) at low temperatures. Also, taking e 
to be the B coupling to the electric current, we should 
replace e — > etan^u,. Due to the fact that the two Feyn- 
man diagrams nearly cancel, the cross-section is highly 
suppressed. A straightforward estimate of the total scat- 
tering rate and comparison to the Hubble rate indicates 
that e would need to be of order 10 5 for \ — v scattering 
to equilibrate at temperatures below m e \ 



(n e <TVrcl) 



T(6Mi 3 ,T) 



(BIO) 



In n e , we include both spin states of electrons and 
positrons. 

Therefore, the ratio 3H/T considered as a function of 
T inverts to give the DM/SM decoupling temperature as 
a function of e 2 M x / This decoupling temperature is 
shown in figure [TT] for several values of <5Mi3. Since we 
consider temperatures near the electron mass, we calcu- 
late the effective species number numerically. This 
includes heating of photons due to annihilation and 
also the neutrino density. 

Since the kinetic mixing is between the B a ^ bosons and 
SM hypercharge and therefore includes mixing with Z^, 
we can ask if the DM comes into equilibrium with neu- 
trinos through the weak force. Above the electron mass, 
scattering from electrons will always dominate scattering 
from neutrinos because mz 3> [i 2 (and because the 
coupling to the weak current is suppressed by ~ A.l m z) 
[35l | . However, below the electron mass, there are many 
more neutrinos than electrons, so neutrino scattering is 
potentially important. Scattering from neutrinos pro- 
gresses through two Feynman diagrams similar to figure 
1101 one with a B propagator and one with a Z prop- 
agator, but there is a relative sign between the two in 
the kinetic mixing. Therefore, to include both diagrams 
properly, we should replace 
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Appendix C: %2 — > Xi e+e a * one loop 



In this appendix, we give a careful derivation of the 
lifetime for \ 2 ~~ ^ Xi e+e ~ decay at one loop level when 
£3 = 0. This is the dominant decay process we have been 
able to find for this case. 



Consider the loops in figure [T|)] As noted above the 
figure, the two diagrams nearly cancel due to the opposite 
signs of the nonabelian 3-gauge-boson couplings; in fact, 
they do not cancel completely only because the two gauge 
bosons in the loop have different masses. The complete 
amplitude can be written as 



M = ig 



Q 2 ~ M§ 



(CI) 



where the spinors of the x particles are labeled by their 
color (as their momenta will be), we take the out- 
going momenta of the positron and electron to be p± 
respectively, and q = k 2 — k\ = p + + p_ . 

The momenta running counter-clockwise around the 
loops are I + 5p on the electron line, /-pon the upper 
gauge boson line, and l+p on the lower gauge line, where 
p = q/2 and Sp = (p + — p_)/2. With these conventions, 
the loop integrals are 
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Vnvfip - l)x + r) v \{2l)n - rj^l + ip)^ 



As usual, we can rewrite the denominators with Feynman parameters as 

dxdydz 2<5(1 — x — y — z) 
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where I' = I +p(y — x) + Spz and we have used p- dp = 0. From the second form, it is clear that taking p\ o p, 2 is 
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the same as swapping the Feynman parameters x -H- y. must come from shifting I to V . 
Therefore, the only terms that survive taking the differ- 
ence in (|C2| must be antisymmetric in x O y, and these In the end, we find 
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with A(x,y) given as in (IC3[) . Note that this loop ap- 
pears to generate an interaction with one extra deriva- 
tive compared to a magnetic moment operator. In the 
approximation, we have taken the gauge boson masses 
to be much larger than any of the momenta in the de- 
nominator. We also have u e 5pv e — 2m e and u e fv e = 0, 
so the amplitude is finally 



with E± and p± the energy and 3- momentum of the . 
Integrating over the nonrelativistic phase space, we find 



M = 
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u\5i>u2u e v e K 
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where A is the function of fj,x,2 in square brackets in (|C4I) . 

In the nonrelativistic limit, the spin-summed and av- 
eraged square amplitude is 



\M\ 2 = 



X A 2 E,E. 



1024 a 2 g a 2 t\t 2 m 2 e M 2 
(E+E- —p + -p- — ml) , 



(C6) 



r = ^^44^ (SM 12 - 2m e ) 3 (5M 12 - m e ) A 2 . 

97T /I 3 

(C7) 

Using the same estimated parameters as below equation 
(f5"3"]) . we find a lifetime of more than 10 33 s. 
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